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ABSTRACT 

We use a hybrid approach that combines high-resolution simulations of the formation of a 
Milky Way-like halo with a semi-analytic model of galaxy formation to study the mass content 
of dwarf galaxies in the concordance ACDM cosmology. We find that the mass within 600 pc 
of dark matter haloes hosting luminous satellites has a median value of ~ 3.2 x 10 7 Mq with 
very little object-to-object scatter. In contrast, the present day total luminosities of the model 
satellites span nearly five orders of magnitude. These findings are in very good agreement with 
the results recently reported in the literature for the dwarf spheroidal galaxies of the Milky 
Way. In our model, dwarf irregular galaxies like the Small Magellanic Cloud, are predicted to 
have similar or slightly larger dark matter mass within 600 pc. 
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1 INTRODUCTION 

The dwarf satellite galaxies of the Milky Way are the most dark 
matter dominated systems known to date in the Universe. They 
represent a heterogeneous population in terms of their stellar prop- 
erties such as luminosity, star formation and chemical enrichmen t 
histories jMated 1 19981 ; iDolphin et all 120051 ; iMartin et alj l2008h . 
Yet, the mass enclosed within a radius of 300 (or 600) pc appears 
to be roughly constant jGilmore et"ai1l2007l : [strigari et alj|2007l 



iThoul & Weinberg!! 19961 : lGnedinll200d : Ibkamoto et alj|2008l, and 
references therein). This has often been considered as a pos- 
sible solution to the excess of small scale structures found in 
CDM an d in particular in TV-body simulations of galaxy size 
systems jKauffmann et aLlll993l; iKlvpin et alj|l999l; iMoore et all 



2008). This could imply a minimum mass scale for the ex istence 
of dw arf spheroidal galaxies, as originally suggested by Mateo 
dl998l) . It is currently unclear whether this is due to the micro- 
physics of the dark matter particles, or to astrophysical processes 
that inhibit star formation on small scales. 

For Weakly Interacting Massive Particles (WIMPs), colli- 
sional damping and free streaming are expect ed to cut off the 
power spectrum at masses of ~ 10 -6 Mq (e.g. lGreen et al.|[2005t) 
or smaller. Then, in models which assume these as dark mat- 
ter particles [e.g. warm dark matter (WDM), A cold dark matter 
(CDM)], a minimum mass scale for dwarf spheroidals can only 
result as a consequence of astrophysical processes that affect the 
collapse of baryons and the formation of stars on small galac- 
tic scales. For example, the presence of a strong photo-ionizing 
background (possibly associated to the reionization of the Uni- 
verse) can suppress accretion and cooling in low-mass haloes. 
This is because the heating of the gas will raise its pressure and 
therefore may suppress its c ollapse in haloes with virial veloci- 
ties < 30 - 50kms _1 (e.g. lEfstathioulll992l : lOuinn et all 1 19961 : 
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ll999] : lBuUocketal.l200ri Somervill e 2002: lBenson et al]2002h . In 
addition, in systems with virial temperature below 10 4 K gas can- 
not cool via hydrogen line emission, and must rely on the highly 
inefficient cooling through collisional exc itations of H2 molecules 
dHaiman et aifeood : iKravtsov et ai]|2004h . 

In this Letter we discuss the existence of a common mass scale 
for Milky Way satellites using results from high resolution TV-body 
simulations of galaxy-size haloes coupled with semi-analytic tech- 
niques to model the evolution of the baryonic component of galax- 
ies. This approach allows us to identify the dark matter substruc- 
tures that host stars and to characterise their stellar properties. At 
the same time, the high resolution of the simulations used in this 
study permits a reliable determination of the internal dynamical 
properties of these satellites. As we shall describe below, we find 
that the dark matter mass within 600 pc for the model satellites 
shows very little scatter from object to object. This is in very good 
agree ment with the observational results by IStrigari et alj [2007, 
2008). Interestingly, our model also reproduces the very wide range 
of luminosities observed for the satellite galaxies around the Milky 
Way. 



2 THE SIMULATION AND THE GALAXY FORMATION 
MODEL 

In this study, we use a high-resolution resi mulation of a ' Milky 
Way' halo from the GA series described in IStoehr et alj j2002t) 
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and IStoehr et alj <2003h . The candidate 'Milky Way' halo was se- 
lected as a relatively isolated halo with a 'quiet' merging history 
(last major merger at z > 2) and with maximum rotational veloc- 
ity close to ~ 220 km s _1 . The halo, selected from an interme- 
diate resolution cosmological simulation, was then re-simulated at 
four progressively higher resolut ions using the 'zoom' technique 
(Tormen. Bouch et & White! Il997l) . The underlying cosmological 
model is a flat A-dominated CDM Universe with cosmological pa- 
rameters: Q m = 0.3, Qa = 0.7, H = 70kms _1 Mpc" 1 , n = 1, 
and as = 0.9. In this study, we use the highest resolution simula- 
tion of the series - GA3new - which contains ~ 10 7 particles within 
the virial radius. 

As explained in |Pe Lucia & Helmil {2008), the simulated halo 
is more massive (Af20o(-z = 0) ~ 3 x 10 12 Mq, - M200 is here de- 
fined as the mass of a spherical region with interior average density 
200 times the critical density of the U niverse at redsh i ft z) than re - 
cent ly estimated for o ur Galaxy, e.g. Battaglia et al. (2005, 2006) 
and IXue et al.1 d2008t> . Q Following iHelmi et alj fe003l) . we then 
scale our 'Milky Way' halo by adopting a scaling factor in mass 
M200/MMW = 7 3 = 2.86. This implies that we scale down the 
positions and velocities by a factor 7 = 1.42. The Plummer equiv- 
alent softening length for the scaled simulation is 0.18 kpc. The 
scaled particle mass is 1.03 x 1O 5 M . 

Simulation data (stored in 108 outputs between z — 37.6 and 
z — 0) were analysed using a standard friends-of-fri ends algorithm 
and the substructure finder algorithm SUBFIND dSpringel et al.1 
l200lh . As in previous work, we have considered all substructures 
retaining at least 20 self-bound particles - which sets the substruc- 
ture mass limit to 2.06 x 10 6 Mq, for the scaled simulation. Fi- 
nally, these halo catalogues have been used to construct the halo 
merger trees that represent the basic input for the galaxy forma- 
tion model used in this study. For de tails on the post-pro cessing, 
we refer the interested reader to ISpringel et all d2005h and to 
|Pe Lucia & Blaizotl fe007h . 

The galaxy formation model used in our study is a refinement 
of the model described in |Pe Lucia & Helmi (2008), who have 
studied the predicted physical properties of the Milky Way and of 
its stellar halo using the same set of re-simulations used here. This 
model buil ds upon the me t hodol ogy introduced in Spring eTet al.1 
d200 lh and |Pe Lucia e t al. (2004) and has b een further refined i n 
later years. The interested reade r is referred to Croton et al] j2006t) . 
|Pe Lucia & Blaizotl J2007I) . and lDe Lucia & Helmil d2008l) for a de- 
tailed account of the modelling of the various physical processes 
considered. In order to obtain a better agreement with the observed 
properties of the Milk y Way satellites, a few refi nements were made 
to the model used in |Pe Lucia & Helmil §008). For completeness, 
we give here a short description of these refinements that will be 
described in detail in a forthcoming paper (Li et al., in preparation). 

• As described in ICroton et al.1 {2006), reionization was mod- 
elled using the formulation provided bv lKravtsov et alj J2004h and 
was previously assumed to start at redshift 8 and to be completed 
by redshift 7. In the study presented here, we assume an 'early' 
reion ization which starts at redshift 15, and is complete by redshift 
11.5 dSpergel et al]|2007l) . 

• In previous studies, haloes with virial temperature below 
10 4 K were allowed to cool at the rates corresponding to 10 4 K. 



1 We note that lBattaglia et alj j2005t) and lXue et a?] <2008l) use a different 
definition of virial mass which is 100 times of the critical density. Our halo 
mass is still larger than observational estimates when the same definition is 
adopted. 
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Figure 1. The solid histogram shows the present-day mass function for the 
satellites in our model. Error bars show the corresponding 1-cr Poisson un- 
certainties. The dashed histogram shows the subhalo mass function, which 
steeply rises up to the resolution limit of our simulation. 



In this stud y, we completely s uppress cooling in these low-mass 
haloes (e.g lHaiman et al . 2000). 

• In previous implementations of our galaxy formation model, 
all metals produced by new stars were instantaneously mi xed to the 
cold gas dDe Lucia et aT1l2004r . Be Lucia & Helmj||2008h. I nspired 
by the numerical simulations of Mac Low & Ferraraldl9991) . in this 
study we assume that for galaxies embedded in haloes with virial 



mass below 5 x 10 Mq, most of the new metals (95 per cent) are 
ejected directly into the hot gas phase. 

As we will show in Li et al. (in preparation), these refinements 
result in a satellite population with physical properties that closely 
resemble those observed for the Milky Way sa tellites, while leav- 
ing es sentially unaltered the results discussed in lDe Lucia & Helmil 
d2008h . 



3 RESULTS 

Fig.Q]shows the mass function of all substructures identified at red- 
shift zero within 280 kpc from the central galaxy and lying in the 
same FOF group (dashed histogram), and the corresponding mass 
function of satellites (i.e. subhaloes hosting stars) in our model 
(solid histogram). The subhalo masses plotted in Fig. Q] are com- 
puted summing up the masses of all self-bound particles, and scaled 
as discussed in Sec. [2] Fig.Q]shows that the dark matter masses of 
our satellites span a relatively large range from ~ 5 x 10 6 Mq to 
~ 5 x 10 10 Mq, and that the distribution is nearly flat between 10 7 
and 10 9 Mq. In contrast, the subhalo mass function continues to 
rise steeply up to the resolution limit of our simulation. The sim- 
ulation used in this study contains almost 2000 subhaloes in the 
considered region, but our model predicts that only 5 1 of them host 
stars. This is still larger (by a factor of about 2) than the number of 
Milky Way satellites currently known. However, when corrections 
due to incompleteness ar e considered, this discr epancy in number is 
eliminated. For example. IKoposov et al.1 J2OO8I) predict ~ 45 satel- 
lites down to My = —5 magn itude. This is consist ent with the 
number of satellites predicted by Tollerud et al. (2008) down to the 
same magnitude limit. The small number of faint satellites could be 
a feature of this particular model, or due to the numerical resolution 
of the simulation adopted in our work. This will be tested in future 
studies using higher resolution simulations. 
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The high resolution of the simulation used in this study allows 
us to measure the satellites' da rk matter mass encl osed within a 
(small) given radius. Following Strie ari et alj d2007t) . we measure 
the mass within 600 pc, which corresponds to 3.33 (scaled) soften- 
ing lengths. We do not attempt to mea sure the mass within 300 pc, 
as done in the more recent analysis by Stri gari et alj {2008), as this 
would be beyond reach for the GA3new simulation. 

In order to measure the mass within 600 pc (Mo. 6 hereafter), 
we compute the centre of mass position for each subhalo using its 
10 per cent most bound particles. We then count the number of 
bound particles located within 600 pc from the centre of mass of 
each subhalo, and multiply this number by the particle mass. We 
find that our subhaloes have on average ~ 380 particles within this 
distance, and in nine cases < 100 particles. 

Numerical effects will in general tend to artificially lower the 
mass in the inner regions of subhaloes. We believe, however, that 
our measurements are robust and are not strongly affected by the 
limited numerical resolution. In order to quantify the impact of nu- 
merical resolution on our estimates of Mo. 6, we assume that the 
inner density profiles of sub haloes are well fit by Einasto profiles, 
as recently demonstrated by ISpringel et al. I J2008h using the very 
high-resolution simulations of the Aquarius project. The Einasto 
profile can be fully characterised with 3 parameters, namely the 
logarithmic slope a, the peak circular v elocity Knax an d the radius 
r mal where this peak value is reached. Springe fet al.1 d2008h find 
q ~ 0.15 — 0.25, with an average value of ~ 0.18. Assuming 
this value for a, and measuring V max and r max directly from the 
simulation used in this study, we obtain an independent estimate of 
M . 6 . 

When we compare the Afo.e values derived using the Einasto 
profile with those directly measured in the simulation, we find at 
most a factor two increase, but in most cases the difference is 
less than 40 per cent. The largest deviation is found for a subhalo 
nearing disruption and whose circular velocity curve is particularly 
noisy. For this object, as well as for those subhaloes with rotation 
curves within 600 pc (and Mo.e) that have evolved significantly in 
the last 2 Gyrs (these are in general the least massive objects), we 
always show both the direct measurement of Mo. 6 and the value 
derived using the Einasto profile. For the other satellites, we keep 
the direct measurements only. It is important to note that the cor- 
rections derived from the Einasto model are smaller than the scatter 
found in the values of Afrj.6 measured for model satellites. 

Fig. [2] shows the measured Mo. 6 for satellites as a function 
of the total F-band luminosity predicted by our galaxy formation 
model. Black symbols show the values of Mo. 6 measured directly 
from the simulations. As discussed above, for satellites with signs 
of tidal disruption we also show the corresponding values estimated 
assuming an Einasto profile with a = 0.18. These are the grey 
symbols linked to the direct measurement by dotted lines. In this 
Figure we have excluded two objects whose stellar masses, as de- 
rived from the galaxy formation model, are larger than the dark 
matter mass of the associated subhalo. These objects are nearly 
fully disrupted in the numerical simulation, but since we do not fol- 
low the tidal stripping of the stars in our model, we could be over- 
estimating their current luminosity by an unknown factor. How- 
ever, we find that only 10 — 15 satellites show evolution in Mo. 6 
and the rotation curves in the last 2 Gyr. We therefore suspect that 
tidal stripping does not affect significantly the luminosity of model 
galaxies. 

Our model satellites span more than four orders of magni- 
tude in luminosity, which is comparable to what is observed for the 
dwarf galaxies around the Milky Way when one includes the re- 
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Figure 2. Mass within 600 pc as a function of the V-band luminosity for 
our model satellites (black and grey solid circles), and for the eight brightest 
dwarf spheroidal galaxies of the Milky Way (asterisks), where the Sagittar- 
ius dwarf has been excluded because it is clearly not in dynamical equilib- 
rium. Black solid circles correspond to Mo. 6 measured via direct particle 
counting. For the model satellites which show clear signs of tidal disruption 
of the associated subhaloes, we also plot the Einasto Mo. a estimates shown 
with grey symbols. 



cently discovered 'ultra-faint satellites'. In contrast, their dark mat- 
ter masses within 600 pc do not differ by more tha n one decade. 
The a sterisks in Fig. [2] show the estimates given by Strig ari et al.1 
d2007h for the classical dwarf spheroidal satellites of the Milky 
Way. Our results are therefore entirely consistent with previous 
analyses claiming the existence of a minimum mass scale of the 
order of 1O 7 M0, and no (or very fe w) subhaloes (in equilib- 
rium) hosting stars below this thr eshold jMateoll9 98; Striga riet al.1 
l2007ll2008l:lGilmore et alj|2007l) . We note that the correlation be- 
tween Afo.e and the V-band luminosity in the model is somewhat 
stronger than in the real data. In the context of our model, satellites 
with smaller Mo. 6 values could exist, but these are expected to be 
strongly tidally disturbed, i.e. out of dynamical equilibrium. 

All our satellites were massive enough prior to accretion to be 
above the cooling limit set by the atomic hydrogen (e.g. ~ 10 9 Mq 
at 2 ~ 1). Objects of such mass should have lost more than 99 
per cent of their mass in order to have Mo.e smaller than 10 7 Mq 
at present time. Note that in Fig. [2] we have also included model 
satellites that could be the counterparts of systems like the Small 
Magellanic Cloud or NGC6822 (i.e. gas-rich). In our model, such 
luminous objects are expected to be embedded in the most massive 
subhaloes at the present time. Hence Fig. [2] shows that we predict 
that the Small Magellanic Cloud should have Mo. 6 only slightly 
larger, namely ~ 10 8 Mq, than the corresponding values measured 
for the dwarf spheroidals surrounding the Milky Way. 

An intriguing question is whether this result implies that the 
Milky Way satellites reside in dark matter haloes of the same total 
mass. We address this question in Fig. [3] where we show the esti- 
mated Mo. 6 as a function of the total dark matter mass given by 
SUBFIND for the objects shown in Fig. [2] While for most model 
satellites Mo. 6 varies in the range ~ 10 7 — 10 8 Mq, their present- 
day total mass spans three orders of magnitude (from ~ 10 7 to 
~ 10 10 M ). Therefore, in our model, a common mass within 
600 pc does not imply the same total mass. This is because this 
scale is generally too small compared to the full extent of dark 
matter subhaloes hosting luminous satellites. In other words, satel- 
lites are not embedded in dark matter subhaloes that can be charac- 
terised by a single parameter, since the relationship between con- 
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Figure 3. Mass within 600 pc as a function of the total mass estimated by 
SUBFIND. The two dashed lines correspond to the upper and lower limits 
for the observational estimat es of Mp .6 for the 8 classical Milky Way dwarf 
spheroidal galaxies (dSphs; Strigari et al. 2007). Symbols are the same as 
in Fig. [I] 



centration and virial mass is not very tight for these objects before 
the time of accretion. 

Our model satellites are found to be embedded in dark mat- 
ter haloes whose p resent total mass is larger than ~ ICTMq. 
Striga rTet alJ J2008h have suggested that this could imply that these 
objects had a characteristic mass of ~ 10 9 Mq at the time of ac- 
cretion (i.e. when they became satellites of the Milky Way). In 
our model, we find a broad range of total dark matter mass at 
the time of accretion for our luminous satellites, that extends from 
~ 3 x 10 8 M s up to 6 x 10 10 M Q . Comparison to Fig.[3]indicates 
that most of the satellites have lost (a significant amount of) dark 
matter due to tidal stripping after accretion. 

The existence of a minimum mass scale for model satellites 
is essentially the result of a combination of the two afore men- 
tioned physical processes: cooling is strongly inhibited for haloes 
with Tvi r < 10 4 K, and reionization prevents the further collapse 
of gas onto low-mass haloes. 

The wide range of total dark matter masses observed today 
reflects in part the initial broad range of masses of the satellites' 
haloes. Furthermore, the tidal field of the Galaxy halo will act to 
increase the difference in present-day mass of these objects de- 
pending on when they were accreted as well as on their orbits 
( iDe Lucia et al.ll2004l ; |Gao et alj|2004l) . For example, our faintest 
objects, which are typically also the least massive today, were ac- 
creted earlier and have suffered significant stripping. On the other 
hand the brightest satellites reside in the most massive haloes, have 
typically been accreted only recently and therefore have not been 
significantly affected by tides. 

A simple argument that illustrates how the measured range 
of halo masses at accretion translates into a wide range of satel- 
lite luminosities is as follows. We assume that the mass accretion 
histories of haloes hosting luminous sa tellites can be descri bed as 
M(z) = M e" 2acZ , prior to accretion dWechsler et~aT1l2002h . M 
denotes the maximum total mass of a halo that is able to grow 
(without being accreted onto a larger structure) until z = 0. This 
expression has been shown to be valid for clusters and galaxy-sized 
haloes but we have checked that it is valid also for the haloes host- 
ing our satellites, with a c ~ 0.2. We also assume that, at the time 
of accretion, each galaxy's lu minosity is propor tional to the mass 
of its dark matter halo (see e.g. IWang et al.l200r3) . If we neglect the 



star formation occurring after accretion, we can write L oc M max , 
where M ma x = Moe _2acZacc . The range of observed luminosities 
at present day can then be estimated as: 

Liow/iu P = [At/ACJ = [M l r/M^\ x e~ 2a ^\ (1) 

In this expression, the first term would be related to the range of 
total masses of the satellites at fixed accretion epoch, and in the 
second term Az denotes the range of accretion redshifts. In our 
simulation, we find M{j ow /M£ p ~ 70 and Az ~ 6, which give a 
range of luminosities ~ 10 3 , not far from what is observed. 



4 CONCLUSIONS 

We have used a high-resolution simulation of the formation of a 
Milky Way-like halo in combination with a semi-analytic model of 
galaxy formation in order to study the mass content of subhaloes 
hosting luminous satellites. The galaxy formation model used in 
this study shows considerable agreement with a large number of 
obser ved galaxy properties (see discussion in IDe Lucia &Helmj| 
2008 and references therein). With the few refinements discussed 
in Sec[2] the same model is able to reproduce quite well the physical 
properties of the observed Milky Way satellites (Li et al. in prepa- 
ration) while leaving essentially unaltered the level of agreement 
with observational data shown in previous work. 

The key result of this Letter is that our model predicts natu- 
rally a common dark matter mass scale within 600 pc for the lumi- 
nous satellites. Our model satellites span nearly five orders of mag- 
nitude in luminosity, while their dark matter masses within 600 pc 
vary only by one decade (between ~ 10 7 and ~ 10 s Mq), in very 
good agreement with recent observational measurements. The to- 
tal dark matter masses of our luminous satellites, however, span 
about three orders of magnitude with the scatter reflecting the lack 
of a tight concentration- virial mass relation, the different accretion 
times, and the different amounts of tidal stripping suffered by the 
parent substructures once they have fallen onto the Milky Way halo. 
The existence of such a scale in the context of our model results 
from the strong suppression of accretion and cooling of gas by low- 
mass haloes after reionization, as well as from the atomic hydrogen 
cooling threshold at T v i r = 10 4 K. These physical processes then 
inhibit the formation of stars in objects that never reached virial 
velocities above ~ 17 km s _ . 

We note that our analysis has been carried out considering the 
mass within 600 pc, and an even tighter distribution may be ex- 
pected when me asuring the mass w ithin 300 pc as used in the most 
recent study bv lStrigari et alj J2008h - This mass scale is beyond the 
resolution limit of the simulation used in this Letter but is within 
reach of the new generation of ultra-high resolution simulations 
like those carried out within the Aquariu s project. 

After the submission of this Letter, iMaccio et al.1 fc009h pre- 
sented a similar study for the satellites of thre e Milky Way-like 
haloes that reproduces the IStrigari et all J2008l) results. Although 
these authors claim a relatively narrow distribution of circular ve- 
locities at the time of accretion, their masses span a similar range 
as ours. They also reach a similar conclusion as we do. Namely, 
that this relation is due to a variety of masses at the time of ac- 
cretion and a broad d istribution of accretion times. More recently, 
IKoposov et al. (2009) has also shown the common mass scale for 
the Milky Way satellites is expected in their CDM based galaxy 
formation models. It is worthwhile to point out that all these mod- 
els have included the effect of a photoionization background and 
accounted for the inefficient molecular cooling below 10 4 K. 
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